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The disk around the Herbig Ae star HD 169142 was imaged and resolved at 18.8 
and 24.5 /xm using Subaru/COMICS. We interpret the observations using a 2D 
radiative transfer model and find evidence for the presence of a large gap. The 
MIR images trace dust that emits at the onset of the strong rise in the spectral 
energy distribution (SED) at 20 //m, therefore are very sensitive to the location 
and characteristics of the inner wall of the outer disk and its dust. We determine 
the location of the wall to be 23I5 ^ from the star. An extra component of hot 
dust must exist close to the star. We find that a hydrostatic optically thick inner 
disk does not produce enough fiux in the NIR and an optically thin geometrically 
thick component is our solution to fit the SED. Considering the recent findings 
of gaps and holes in a number of Herbig Ae/Be group I disks, we suggest that 
such disk structures may be common in group I sources. Classification as group 
I should be considered a support for classification as a transitional disk, though 
improved imaging surveys are needed to support this speculation. 

Subject headings: circumstellar matter — stars: pre-main sequence — protoplan- 
etary disks 



1. Introduction 



Transitional disks are a class of protoplanetary disks whose inner regions are devoid of 
small dust grains. These disks have been attracting attention recently, since such cleared 
inner regions or gaps may be related to on-going planet formation. The presence of an 
inner hole /g a p has been suggested for the disk around the Herbig Ae star HD 169142 by 
Grady et al.l (120071 ) and iMeeus et al.l (120101 ). mostly based on the analysis of the spectral 
energy distribution (SED). Both groups conclude that the steep rise in fiux at ~20 /im in the 
SED refiects emission from the wall of the outer disk. The location of the wall (i . e. the inner 
edge of the outer disk) h a s not been determined conclusively - iMeeus et al.l (120 lOf ) used 
20 AU while lOradv et al.l (120071 ) suggested 44 AU. Observations of HD 169142 at various 
wavelengths have so far not been able to constrain the inner cavity radius directly and 
accurately. 

One of the limitation of SED modeling is the inability to locate and prove the existence 
of gaps in proto-planetary disks. For instance, gaps may not be reveal ed in the SED be cause 
the inner disk can partially obscure large parts of the outer disk ( lAcke et al.l |2009| ) and 
gaps hidden under that shadow do not leave a fingerprint in the SED. Clever SED modeling 
can reproduce SEDs of truly transitional disks without the need of a gap. The transitional 
disk LkCa 15 has an outer disk starting at 46 AU as seen at 1.4 and 2.8 mm wavelengths 
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( iPietu et al.ll2006[ ). However, llsella et al.l (|2009[ ) show that a model approach with a smooth 



distribution of material from a few stellar radii to about 240 AU can sometimes reproduce 
both the observed spectral energy distribution and the spatially resolved continuum emission 
at millimeter wavelengths for this object. Only from direct imaging it is possible to directly 
constrain the radial density structure of dust in a proto-planetary disk. 

Dust in proto-planetary disks covers a temperature range from ~1500 K at the inner 
dust sublimation radius to a few Kelvin in the outer parts of the disk. Direct mid-infrared 
(MIR) imaging at 18.8 and 24.5 ;um is most sensitive for 100-150 K blackbody dust, though 
strong contributions from dust elsewhere in the disk at lower and higher temp eratures can 



also b e present. If the disk has a very strong NIR component such as HD 135344 flMeeus et al. 



200ll ). the inner region may still be strongly represented in the center of the image at 18 



micron. On the other hand, for disks with a larger gap or an inner hole, the image size 



will b e a sensitive tracer to the location of the inner radius of the outer disk (jVerhoeff et al. 



201l|). 



In this paper, we present direct MIR imaging observations of HD 169142 at 18.8 and 
24.5 /xm, using COM ICS on the 8.2m Subaru telescope. We confirm the result found by 



Marinas et al.l ( 1201 ll ) and find that the disk shows extended emission at 18 /im. In addition 
we find that the disk is also resolved at 24.5 yum. By constructing a radiative transfer disk 
model that fits both the SED and the imaging results, we find that the size of the source at 
MIR wavelengths is most naturally explained by a disk with a large inner gap, i.e. transitional 
disk. We discuss the structure of the disk, and the implications for its nature. 



Observations and Data Reduction 



Obse rvations were condu c ted using COMICS (Cooled Mid- Infrar ed Camera and Spec- 



trometer; iKataza et al.l l2000l : lOkamoto et al.l l2003t ISako et al.l |2003|) on the 8.2 m Sub- 



aru Telescope on Mauna Kea, Hawaii. HD 169142 was observed using the Q24.5-NEW 
(A=24.5 /xm, AA=0.75 /im) and Q18.8(A=18.8 /im,AA=0.9 fim) filters. The plate scale of the 
COMICS camera is 0.13 arcsec per pixel. The chopping throw was 10 arcsec and the posi- 
tion angle of chopping direction was degree. The chopping frequency was 0.45 Hz. The 
total integration time for Q24.5 and Q18.8 observations were 802 s and 360 s, respectively. 
Just before and after observing HD 169142, we took data of PSF reference and photometric 
standard stars. We used 6 Oph for Q24.5 and a Her for Q18.8. The total integration times 
for the reference stars were 243 s and 83 s, respectively. A summary of the observations is 
given in Table 1. 
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We processed the data using a shift-and-add method to improve the blurring caused 
by atmosphere, tracking errors, misregistration, etc. The imaging data consists of 0.98 s 
on-source integration frames. First, the thermal background and the dark current signals 
were removed through the subtraction of the chopped pair frames. The object and the 
PSF star are bright enough to be recognized even in 0.98s chop- subtracted frames, so we 
searched the centroid of the object. We then shifted the frames so as to align the centroid 
position, and summed up the frames. However, we excluded frames which were blurred by 
atmospheric seeing. We rejected frames whose radius of 78.4% encircled energy (r784) is 
larger than thrshould radius (rc) which covers 95% of the wavefronts whose Strehl ratio is 
0.9. These thrshould values are determined as = 1.131" at 18.8 ^m and = 1.482" at 
24.5 yum by Monte Carlo simulation of ideal unresolved point source. Due to this rejection 
of lower quality data, the effective integration time of HD 169142 was reduced to 345 s and 
215 s for Q24.5 and Q18.8, respectively. The same procedure was applied to the PSF images 
taken before and after HD 169142, where the effective integration time became to 195 s and 
83 s, respectively. Since HD 169142 shows marginal extention compared to the PSF and the 
rejection criterion Tc is determined for the ideal unresolved point source, frame rejection 
rates of HD 169142 are higher than those of PSF stars. It means that the effective rejection 
criterion is slightly more stringent for HD 169142 than for PSF stars, however, it does not 
overestimate the extention of HD 169142. 



For flux calibration we used template spectrum provided by lEngelke et al.l ( 120061 ). Using 
the standard stars observed over the two nights, an airmass correction was applied for Q24.5 
photometry. For Q18.8 photometry, we could not find a significant airmass dependence 
during the night. A standard aperture photometry was applied and the resultant flux density 
of HD 169142 was 10.5 ± 0.4 Jy at 18.8 /xm and 13.0 ± 0.5 Jy at 24.5 fim. The final images of 
the PSF and HD 169142 are shown in Figure [H 



3. Results 

3.1. Source size of HD 169142 in 18.8 and 24.5 fiin 

The azimuthally-averaged radial brightness profiles of HD 169142 and the PSF stars at 
18.8 and 24.5 are shown in Figure [2l It is clear that HD 169142 is extended at these 
wavelengths. At 18.8 fim the direct FWHM of the shift-and-added images of HD 169142 
and of the PSF stars were 0.604" ±0.017" and 0.493" ±0.006" , respectively. At 24.5 ^m, the 
corresponding numbers are 0.680" ±0.034" and 0.628" ±0.007" (see Table 1). It is surprising 
that the size of HD169142 is not increasing with wavelength. For a cont inuous flarin g disk. 



one would expect the size of the image to scale roughly with the PSF f lMeijerl 120071 ). Our 
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images therefore give a first indication that the radial density structure is not continuous. 
The FWHMs of the PSF references are comparable to the predicted value of the diffraction 
limited performance of the telescope. As a rough estimate of the intrinsic source size of 



HP 16 9142, we applied the quadratic subtraction method as described by iMarinas et al. 



( 1201 ih . The de rived source F ^ yHM s were 0.349" ±0.014" at 18.8 //m and 0.261" ±0.025" at 
24.5 /im. Since iMariiias et all J201lh observed the FWHM of this source to be 0.32" ±0.05" 
at 18.0 /im, our measurement at 18.8 /xm is in good agreement with their results within the 
uncertainties. 



4. Modeling 



4.1. Observational constraints 



Complementary views of the disk structure and the dust properties of HD 1 69142 are ob 



tained by adopting photometric data at various wavelengths from the 


iterature 


van den Ancker et al. 


1997; 


Cutri et al. 


2003; 


Zacharias et a .. 


2004; 


Svlvester et al. 


1996; 


Meeus et 


al. 


2010 


) and 



MIR spectra from ISO (IMeeus et al.l l200ll ) and Spitzer/IRS (jJuhasz et al.ll2010l ). These 
datasets reflect the disk brightness as a function of wavelength and are a result of the disk 
structure defined by the density distribution, composition, inclination, and inner and outer 
radius of the disk and by the pr operties of the central star. An outer disk radius of 235 AU 
is taken from lPanic et al.l (120081 ). An inner disk radius of 0.1 AU is set at the l ocation where 



the d ust reaches 1500 K and starts to sublimate. With an inclination of 13° (IRaman et al. 



20061 ). the system is nearly pole-on and, therefore, deviations from axial symmetry caused 
by projection effects can be neglected in our images. Th e stellar spectrum is described by a 
Kurucz model with an effective t emperature of 8200 K (IDunkin et al.l 119971 ). luminosity of 
15.3 LfT^ and extinction At/=0 . 46 (Ivan den Anckerlll999l ). The system is set at a distance of 



145 pc (jde Zeeuw et al.lll999l ). All stellar and disk parameters used in this paper are shown 
in table H 



4.2. Radiative transfer code MCMax 



We construct a geometrical model of the disk of HD 169142 and fit it to both the 
radial surface brightness profile (RBP) of our obs erved images an d to the SED. For the 
model we use the radiative transfer code MCMax (IMin et al.ll2009l ). which is able to solve 
temperature and density structures in very opticall y thick circurnstella r disks. This code has 
been compared to other radiative transfer codes by iPinte et al.l (120091 ) and has been applied 
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successfully in previous studies (e.g. IVerhoeff et al.ll201ll : iMulders et al.ll201ll ). We assume 



an axisymmetric dust distribution in which t he optical proper ties of the grains are computed 
using a distribution of hollow spheres (DHS; iMin et al.ll2005[ ). We consider all dust species 
in the disk to be in thermal contact, and calculate their respective opacities from the optical 
constants as if they were separate particles. Furthermore we assume the gas temperature 
to be set by the dust temperature. We adopted grains with sizes acco rding to the power- 
law /(a) oc a~^'^. The power-law index is that of interstellar grains (jMathis et al.l 119771 ) 



and is consistent with collisional fragmentation considerations ( iHellyerl 1 19 701 ). Under the 
assumption that the system is in hydrostatic equilibrium, the vertical density distribution 
of the disk is solve d by iterating the densi ty and temperature structures until they become 



self-consistent (e.g. iDuUemond et al.l 120071 ). 



We assume a gr ain population consisting of 30% carbon and 70% silicates (we refer to 



Mulders et al.l (120111 ) for references of the optical properties of this dust composition). The 



contribution of warm small (< 10 fim) amorphous silicate grains to the spectrum must be 
low, since the flux level at 10 ^m is very low compared to that at 20 /xm. We model this 
change in flux level by replacing the amorphous silicate dust in the inner disk with dust with 
a higher continuum opacity. We here use carbon, but other d ust species such as met allic 
iron may give the same result. This treatment is consistent with Ivan Boekel et al.ll2004l who 
find that grain growth in the innermost regions has proceeded further than in the outer disk 
regions. We have not attempted to fit the PAH features in detail, though we implement 
them into our model to fit the SED. Scattered light contribution is also included and the 
fraction of the scattered light to the total flux density in 18.8 and 24.5 /im is 22-26% in this 
model. 



4.3. The best fit model 

We fit the SED as well as the 18.8 and 24.5 yum images simultaneously by using the disk 
parameters shown in Table |2i The RBPs and SED of our final best fit model are shown by 
the the purple dashed lines in Figures [2] and [31 A cartoon sketch is given in Figure |H In 
this section we will first discuss the geometrical characteristics of the location of the gap and 
the outer disk properties. Thereafter we discuss the constraints we find to the inner disk 
structure of the best fitting model. 
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Fig. 1.— Shift-and-added image of HD 169142 (top) and PSF star (bottom) at 18.8 and 
24.5 ^m. Brightness unit is in Jy/arcsec^ and scaled from to the peak value in the image. 
North is up and the east is to the left 
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Fig. 2. — Left and middle panels are model brightness distribution of final model and model 
image after convolved by Subaru/COMICS PSF (middle), respectively, shown in Jy/arcsec^. 
Right panels represents the peak-normalized azimuthally-averaged radial brightness profiles 
relative to the centroid of the image (right). Top panels are in the Q18.8 and the bottom 
panels are in the Q24.5 filters. In the radial profile plots (right), the observation of HD 169142 
is shown by the red diamonds, and the corresponding PSF of the calibrators is shown by the 
black squares. The final model convolved with the PSF is shown by the black dotted line. 
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Fig. 3. — Spectral energy distribution of HD 169142. The solid grey lines represent the 
Kurucz and the Spitzer/IRS [5.3 - 37 ^m] spectra respectively. The dark green dashed line 
shows the geometrically thick, optically thin halo-like dust cloud component. The red dotted 
line gives the flux contribution from the innermost part of the outer disk (23-26 AU). The 
flux from the remainder of the outer disk (26-237 AU) is given by the dashed-dotted blue 
line. The sum of the components is given by the purple solid line. It is clear that at both 
18.8 and 24.5 ^m, the SED is dominated by flux coming from the inner edge of the outer 
disk. 



4.3.1. Fitting the SED 



Previous SED modeling by lMeeus et al.l ( l2010l ) and iGrady et al.l (120071 ) suggested a gap 
in the disk with an outer radius of 20 and 44 AU respectively. We confirm and improve 
on this result by our test modeling for HD169142: disks without a gap fail to fit the SED, 
despite various attempts within the physically reasonably parameter space. So already the 
SED fitting gives a str ong indication for the presence of a gap (confirming the results by 
Meeus et al.l (120 lOf ) and lGrady et al.l ( l2007l ). even though its size is not well constrained at all 
from the SED alone. We incorporate this result and include a gap in our model. The effect 
of the discontinuity in the radial density structure is that a large fraction of the stellar light 
now falls on the vertical wall at the inner edge of the outer disk. Emission from this wall is 
clearly distinguishable in the SED (Figure [31 illustrated by the red dotted-line). In addition. 
Figure [3] shows that the flux contribution from the inner edge (between 23-26 AU) dominates 
at 18.8 and 24.5 /im and that flux contributions from the inner regions and from the outer 



- 10 - 




J_J 1 I I , I I 1 I 1_J 

0.1 1.0 10.0 100.0 

AU 



Fig. 4. — Sketch of our final model with two possible structures for the inner regions. Top: 
the model with a flat inner disk, a optically thin halo and an outer disk. Bottom: the model 
with a denser but still optically thin halo, and an outer disk in which the scaleheight of the 
dust is decreased. Both geometries give an equally good fit to the SED and the images. 
Although we can not distinguish between these models with our dataset, we choose the 
bottom model as our final model described in the text. 

disk are negligible at these wavelengths. Therefore, our Subaru/COMICS images at 18.8 and 
24.5 fim are sensitive tracers of the radius where the outer disk starts. Another insight to 
understand the dominance of the wall is given by the radial density and temperature plots of 
the outer disk of this model in Figure 13 The dotted lines show that the densest regions are 
largely optically thick. Therefore the temperature of the disk mid-plane decreases rapidly as 
a function of radius. Thus 18.8 and 24.5 fim photons are mostly produced at the inner few 
AU radii of the outer disk. 

The radius where the outer disk starts is not well constrained by SED modeling. Be- 
cause the disk structure is degenerate by uncertainties of many parameters like the stellar 
temperature and luminosity, the interstellar extinction, the grain size population of the dust, 
the vertical density structure and so on. Only by fitting our model to the resolved MIR ob- 
servations we can unambiguously confirm the existence of a gap and constrain the location 
of the wall of the outer disk. 
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Fig. 5. — Left: the dust density pdust in g cm~^ of the outer disk as a function of radius 
and angle 6 for our best fitting modeL The lower and upper density cut-offs are 10"^" and 
10~^ g cm~^ respectively. Right: the temperature structure of the outer disk as a function 
of radius and angle 6 for our best fitting model. The inner edge of the outer disk starts at 
a temperature of 180 K. For small 6 the angle is approximately Z/r. The grey dot-dashed 
line shows the r=l surface in radial direction in scattered light (A=0.5 /im). The dashed 
lines represent the r=l surfaces as seen perpendicular from the disk. From top to bottom 
the colors are red, purple, orange and blue and represent the wavelengths 0.5, 18.8, 24.5, 70 
/im. 



4-3.2. Fitting the MIR-images, constraining the location of the wall 



After convolving our best fit model with the PSFs of the calibrators, we find that the 
imaging data is fitted best with an inner edge of the outer disk of 23l;5AU (see Fig|2]). 
The uncertainty on this radius is determined by fitting the uncertainty of the azimuthally 
averaged flux level at ±1 sigma. This result is consistent with a wall radius of 20 AU 



( iMeeus et al.ll2010[ ) but excludes the model with a wall at 44 AU (IGrady et al.l 120071 ). 



We have performed a small parameter study where we have examined the robustness 
of our final model to fit the RBPs at MIR wavelengths. We have varied the following 
parameters: the surface density powerlaw between {-0.5, -1.5}, the grain size distribution 
powerlaw {-2.5, -4.5}, the dust settling factor {0.5, 1}, the stellar temperature {8000 K, 8500 
K} and luminosity {13 Lq, 17 Lq} and the distance {130 pc, 160 pc}. While the RPB only 
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vary slightly in this parameter study (all within the 1 sigma error of the measurements), the 
SED fit becomes significantly worse. We conclude that the radius of the wall is the only free 
parameter which significantly influences the broadness of the RBP. Therefore our fit to the 
RBP is robust in constraining the radius of the inner edge of the outer disk at '22^^^^ AU. 



J^..3.3. The inner disk 



In this section we explain in more detail the implications for the inner disk structure. 
Modeling shows that an optically thick hydrostatic inner disk does not produce enough flux 
to fit the SED. Scaling up the vertical height of the inner disk would not be a solution because 
if the inner disk is higher than Z/r ~ 0.1 it covers the outer disk for more than ~50% in 
direct and scattered light. This inner disk then casts a too large shadow on the outer disk to 
be able to reproduce the SED at wavelengths longer than ~10 fj,m. To solve this problem we 
include a geometrically thick, optically thin halo-like dust cloud to the inner regions (< 0.5 
AU, see Figure E]) to fit the SED in the NIR and to avoid a shadowing effect on the outer disk. 
The physical origin is not certain at this moment, though other observational studies show 



that an optically th i n component close to the star provides a good fit as well ( IGrady et al. 
200?! : benistv et aD boili IVerhoeff et aP boili iMulders et all boiol). F urthermore, fitting 



the Herbig stars' median SED requires such a halo ( iMulders fc Dominikll2012[ ). Further high 
spatial resolution observations such as infrared interferometry will be necessary to reveal the 
inner disk structure. 

If there is no shadowing effect from the inner disk, the wall of the outer disk is too 
luminous. To fit the SED again, we can either settle the grains in the outer disk by scaling 
down the dust density scaleheight by a factor of 0.6 or we can scale down the abundance 
of small ( < 1 /um) grains by one o rder of magnitude. These ad justments are consistent 
with dust settling and grain growth ( iDuUemond fc Dominikll2005l ). however, it is uncertain 
to what extent these effects are present in the disk. The influence of other effects like 
dynamical scattering of larger dust grains may also play a role in the vertical grain size 
distribution. Alternatively, we can combine a halo-like dust cloud and a flat inner disk and 
consequently perform no settling or removal of small dust grains. So we conclude that there 
are two possible inner disk geometries (as illustrated in the cartoon in Figure H]) which give 
an equally good flt to the SED and the RBP. 

Without an optically thick inner disk, the RBP is still reproduced with the inner radius 
of the outer disk at 23l^ AU. The flnal model can be interpreted as an inner halo-like dust 



cloud ( iKrijt fc DominikI 1201 if ) and a gas-rich outer disk structure, so truly a transitional 
disk. 
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Discussion 



5.1. Inner hole/gap are common for Meeus's group I sources ? 

We have shown observational evidence for the presence of an inner hole or gap in 
the HD 169142 disk. A similar st ructure is als o indicated for other Herbig Ae/Be stars 



such as AB Aur Jflonda et al.l boioh. HP 142527 JPukagawa et 



Verhoeff et al IbnilL HD 135344 JBrown et al.ll2nn9h. HP 36 112 



2006 



Fuiiwara et al.ll2006 



■I boioh 



f lBouwman et al 



20031 : iBenistv et al 



2010 



Isella et al 



HD 100546 



Mulders et al.l 1201 if ) . All of these objects have 



been classified as group I according to the definition introduced by lMeeus et al.l (120011 ). Since 
the strong far infrared excess is a characteristic of group I sources, we speculate that this 
component may come from the emission from the inner-edge of the outer disk. Consequently, 
we suggest that the inner hole or gap structure may be common for group I sources. If this 
is true, the hypothesis that the group I fiared disks evolves into the group II fiat disks due 
to dust sedimentation should be reconsidered, because the difference of group I and II is not 
only the degree of disk fiaring or grain growth, but also the presence of cleared inner regions 
(holes and gaps) in the disk. 

As we have shown in the above analysis, the spatial extent of a Herbig Ae/Be disk 
at MIR wavelengths, often described by a FWHM, is strongly infiuenced by the position 
and temperature of the wall. Previous studies showed that group I sources are likely to be 
more extended than group II sources in the MIR wavelengths, and suggested that the disk 
geometry (fiaring or flat disk) may play an important role in the thermal stru cture and MIR 
emission of the disk ( Leinert et al.ll2004j : iLiu et al.l 120071 : iMarifias et al.ll201ll ). Our 24.5 /im 
imaging survey of Herbig Ae/Be stars with Subaru/COMICS (Honda et al. in preparation) 
also conflrms that many group I sources are extended, however, their spatial extent shows 
a great diversity, from marginally to remarkably extended. Such diversity in the spatial 
extent in the MIR can be understood in terms of the distance, the inner disk structure and 
the location of the wall of the outer disk. In fact, the temperature of the wall both in the 
AB Aur and HD 142527 systems appears to be relatively cool (70-100 K) and its wall radius 
also tends to be at some distance away from the central star ( 100 AU in AB Aur; 170 
AU in HD 142527 jFujiwara et al.ll2006l )). The wall in HD 169142 is at a smaller distance 
(23AU), thus the spatial extent is not so large compared to that of AB Aur and HD 142527. 
We suggest that not only the geometry (flared or flat) of the disk, but also the inner gap and 
a wall-like inner edge of the outer disk are important for understanding the spatial extent of 
Herbig Ae/Be disks in the MIR. 
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5.2. Origin of the inner hole of HD 169142 disk 



Gaps seem to be a common characteristic of the group I flaring disks ( 


3onda et al. 


2010 




Fukaeawa et al. 


2006: 


Fuiiwara et al. 


2006: 


Verhoeff et al.l 


2011 




Brown et al. 


2009: 


Isella et al. 


2010; 


Bouwman et al. 


2003 


: Isenistv et al. 


2010l 


Mulders et al. 


201] 


.). An 


inner 



hole or gap in the disk has been explained by several mechanisms, in cluding (1) photoevapo- 



DuUemond & Dominik 


zone ( 


van Boekel et al. 



2004), (3) grain growth in the inner disk causing a lower dust opacity zone ( Ivan Boekel et al. 
20051 ). and (4) the presence of another body i n the disk that dynamically cr eates a gap and 
decouples the inner disk from the outer disk ( lAugereau fc Papaloizoull2004l ). We prefer so- 
lution (4) since this also provides a natural explanation for a inner halo-like dust cloud in 
the inner disk. The presence of a halo likely requires some dynamical interactions with 
planetar y bodies since new gen erations of dust can than be created in situ up to large disk 
heights ( iKrijt fc Dominikll201ll ). Grain growth can not explain the presence of a gap and 
a wall. Furthermore, our test models show that geometric shadowing (i.e. varying scale- 
heights, surface density powerlaws) of a part of the disk at several AU can not reproduce 
the SED and the RBP of the images. We recognize that planetary companions have not yet 
been detected in the transitional disk around H D 169142, though rec ent dir ect detection of 



the pl anets around A type stars (e.g. HR8799; iMarois et al.l ( l2008l ). dPic iLagrange et al. 



( I2OIOI )) leads support for such hypothesis. Thus we suggest that the HD 169142 system is 
an excellent candidate to look for newly-formed planets in the protoplanetary disk. 



6. Conclusions 

We arrive at the following conclusions: 

1. HD 169142 is extended at MIR wavelengths with the Subaru/COMICS instrument. 
The derived source sizes from quadratic PSF reference subtraction are 0.349" ±0.014" 
at 18.8 ^m and 0.261" ±0.025" at 24.5 /xm, respectively. 

2. The observed sizes as well as the steep rise of the SED near 20 nm require a disk model 
with an inner hole the size of 23^5 AU, making HD 169142 a transitional disk. 

3. HD 169142 is one of a growing number of Herbig Ae disks that are classified as group I 
source and that also require the presence of a large inner gap. It appears that many if 
not all group I sources may be strong candidates for classification as transitional disks. 
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Table 1. Observations Summary 



object 


filter 


Date(UT) 


Integ. time [s] 


AirMass 


Direct 








used/total (%) 




FWHM 


50ph 


Q24.5 


2004/07/11 


195/243 sec (80 %) 


1.254-1.354 


0.628" ±0.007" 


HD 169142 


Q24.5 


2004/07/11 


345/802 sec (43 %) 


1.541-1.542 


0.680" ±0.034" 


aHcr 


Q18.8 


2004/07/12 


83/83 sec (100 %) 


1.081-1.115 


0.493" ±0.006" 


HD 169142 


Q18.8 


2004/07/12 


215/360 sec (60 %) 


1.553-1.559 


0.604" ±0.017" 



Table 2. Parameters of HD 169142 system used in our best-fit model 



parameter 


value 


Remarks 


Spectral type 


A5Ve 


Dunkin et al. ( 1997) 


Extinction A^ 


0.46±0.05 


van den Ancker fl999) 


log g 


4.22 


van den Ancker fl999") 


Temperature 


8200K 


Dunkin et al. ( 1997) 


Distance 


145±15 pc 


de Zeeuw et al. (1999) 


Age 


6+1 Myr 


Grady et al. ('2007) 


Stellar luminosity 


15.33±2.17 L© 


van den Ancker (1999) 


Stellar Mass 


2.28±0.23 Mq 


van den Ancker (1999) 


Stellar Radius 


1.94±0.14 Ro 


van den Ancker (1999) 


Gas disk mass 


(0.16-3. 0)xl0-2 Mq 


Panic et al. (2008) 


Dust disk mass 


4x10-'' Mq 


fit to the sub mm photometry 


Inclination 


13° 


Raman et al. (2006): Dent et al. (2005) 


accretion rate 


<10-8Mq yr-1 


Gradv et al. (2007) 


Rhalo 


0.1-0.2 AU 


geometrically heigh, optically thin component to fit the NIR 


Rin 


23^5 AU 


fit to RBP of Subaru/COMICS data 


Rout 


235 AU 


Panic et al. (2008) 


surface density exponent 


-1.0 


Hydrostatic equilibrium 


Particle size 


a ={0.03/im,lcm} 


powerlaw distribution of -3.5 


Silicates 


70% 


similar to Mulders et al. (2011) 


Amorphous carbon 


30% 


Zubko et al. (1996) 


MpAH 


0.45 X lO"'^ Mq 


uniform PAH distribution 


Mhalo 


0.28 X 10-1'' Mq 


only carbon 




0.3 X 10-^ Mq 


mass of grains a ={0.03/im, 1cm} in the disk 



